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ASCENT FROM THE LUNAR SURFACE 


By 

Rowland E. Burns 
and 

Larry G. Singleton 


SUMMARY 

The problem of three-dimensional optimal ascent from the lunar surface is discussed in the 
report using the techniques of variational calculus. The Moon is assumed to be spherical and ro- 
tating, but perturbational effects from all other bodies are neglected. Final orbital inclination is 
calculated under the assumption that the angular displacement of the moon is negligible during 
ascent. Only single stage vehicles are considered and are subdivided into propellant and final 
mass in orbit, i.e. , there is no consideration of payload dependencies upon structural mass, etc. 
Furthermore, both the thrust and mass flow rate are assumed constant throughout the powered 
flight for a given vehicle. 


SECTION I. EQUATIONS OF MOTION 

The equations of motion will be derived in a spherical coordinate system. As shown in 
FIG 1, 0 is the latitude angle (measured from the equatorial plane positive toward the north) 

and 0 is the longitude angle (measured from <he Moon's prime meridian positive in the direction 
of rotation). Let r be the radius vector from the center of the Moon to the vehicle. At the tip 
of the radius vector let us define an orthogonal system of unit vectors, u f , uq , u <£. Let u T 
be collinear with r (positive in the same direction as r); let uq be perpendicular to u r in a 
plane containing u T and the polar axis (positive in the direction of increasing Q ): let u $ be 

perpendicular to both « f and uq and chosen to form a right-handed system with u t and uq. 

By trigonometry we can express u t , uq , and u ^ in terms of i , j , and k as 

u T = cos 0 cos (f) i + cos 0 sin <$> / + sin 0 k 

uq - — sin 0 cos <f> i - sin 6 sin <f> j + cos 0 k 

uj, = — sin 0 i + cos <£ / 


Equations (1) through (3) may be differentiated to yield 


( 1 ) 

( 2 ) 

(3) 
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u q = — 6 u r — <f> sin d u ^ 

«-0 = — (cos 0 u T — sin 6 ug ) 

The radius vector, r, may be expressed as 
t = r u r 

The velocity may now be defined as 

. . ^ x . ^ • 

v = r = t u f + r u r = r u r + r 0 ug + r $ cos 6 

The acceleration is given as the time derivative of velocity as 


(5) 

( 6 ) 

(7) 

( 8 ) 


a. = \ = r u r + 2 r u r + r u r = 

[ r — r ( d 2 + cj> 2 cos 2 6) ] u t + [ r 6 +2 rd + r <£ 2 sin 6 cos 0] ug 

+ [ r <j) cos $ — 2 r 0 0 sin 0 + 2 r <jS cos 6 ] u <£, (9) 

The thrust vector, shown in FIG 2, has direction cosines y and 8. From FIG 2 we may 
write the thrust as 

T — T cos 8 cos y u z + T sin 8 u g + T cos 8 sin y u ^ (10) 

The remaining force, gravity, is always directed along the radius. If we denote the accelera- 
tion of gravity as g, the mass of the Moon as M and the universal gravitational constant as G 
we can write 

MG - 

8 = ~~ —yz u * ( 11 ) 

By Newton’s second law we may now write 

m a = T + m g (12) 

where m is the mass of the vehicle. 


Inserting equations (9), (10) and (11) into (12), equating components, and dividing by m 
gives 

t — r ( d 2 + d> 2 cos 2 6) = — cos 8 cos y — ■ (13) 

r m T * 

..... f 

r 6 +2 r 6 + t <f> 2 sin d cos 6 = — sin 8 (14) 

m 

r d> cos 0 - 2r d> $ sin 0 + 2 r d> cos 6 = cos 8 sin y * (15) 

m ' 

We may now state the problem as follows: Subject to constraints (13), (14) and (15), as well 
as specified initial and final conditions, determine the functions y and 8 as functions of time 
such that the maximum payload is injected into orbit. 


* An alternative derivation of these equations is given in reference 4 . 
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The solution of this problem requires the use of the calculus of variations. As equations 
(13), (14) and (15) now stand we would have to make use of second order Lagrange equations. 
To avoid this, we make the following kinematical substitutions: define 


r = p 

(16) 

b = (o 

(17) 

ii 

(18) 

Substituting (16), (17) and (18) into (13), (14) and (15) and isolating all derivative terms 
gives 

P — — cos 8 COS y + ^ G _ r ^2 + a 2 cos 2 Q) _ Q 

m r 1 

(19) 

(6 — — sin 8 + — + a 2 sin 6 cos 6 - 0 

mr r 

(20) 

a — — — -pr cos 8 sin y — 2 co a tan 6 + = 0 

m r cos v r 

(21) 

One further relationship is required. Since the vehicle is assumed to 
flow rate, m, we may write 

have a constant mass 

m — ttiq — m t 

(22) 

where m is the instantaneous mass, mo the initial mass, and t is 

the time since lift-off. 

SECTION II. BOUNDARY CONDITIONS 


Before proceeding to the actual variational formulation of the ascent , it will be advan- 
tageous to determine the boundary conditions associated with the problem. We shall denote 
initial values by a zero subscript and final values by an / subscript. At the initial point, 
t — t 0 = 0, the vehicle will be assumed to be at rest on the lunar surface. The launch longitude, 
0 O , is arbitrary since we have assumed a spherical Moon. The simplest assumption on <^ 0 i s 
that it be chosen equal to zero. Furthermore, the first derivative of cf)Q, i s simply given 

by the angular rate of rotation of the Moon which we shall denote by 12. The latitude of launch, 
(9 0 , is arbitrary.* Since u $ and uq are perpendicular, the lunar rotation induces no motion in 
uq direction and we may set 0 Q = 0. 

In summary, we have the following set of initial conditions: 


t = * 0 = 0 

(23) 


* Except for ± — . The coordinate system chosen is degenerate with respect to 0= ±— . 

& 

If launches into polar orbit are of interest, the rotations which produce equations (1), (2) and 
(3) may be reversed in order of application. The resulting equations will be degenerate with 
respect to equatorial launch, d Q = 0, but acceptable to study launch into polar orbits. 
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(24) 


m = m Q 


r = r 0 = r t (radius of the Moon) 

(25) 

r = r Q = po = 0 

(26) 

6 = 6 0 (arbitrary, but ^ ) 

(27) 

9 = 0 O = co 0 - 0 

(28) 

O 

ll 

0 

-e- 

11 

-e- 

(29) 

= <j> 0 = ct 0 = ft (angular velocity of the Moon) 

(30) 


At the final point we require a circular orbit of specified inclination. From equation (8) we 
may write the square of the velocity as 


v • v = r • r = r 2 +r 2 ( <9 2 + <f> 2 cos 2 6) = p 2 + r 2 {(o 2 + a 


0) 


(31) 


The condition for circularity may be stated as 


r • r = 


M G 


(32) 


Thus 

p 2 + t 2 {co 2 + a 2 cos 2 9) = —— ■ (33) 

r r 

or 

M G — r p 2 — r 3 (n> 2 + < 7 2 cos 2 0) = 0 (34) 


But, for a circular orbit we must also have at cutoff 

; = P = 0 (35) 

Thus (34) becomes 

MG — r 3 ((U 2 +(t 2 cos 2 0) = 0 (36) 

Equations (35) and (36) must be’ fulfilled at / = t f . 

We now digress slightly to determine the orbital inclination as a function of the cutoff para- 
meters. 

Assuming that the rotation of the Moon is negligible during ascent,* the final inclination is 


* The actual amount of rotation is less than two tenths of a second of arc for the trajectories 
considered. Reference 5 presents a detailed treatment of the effect because of rotation of the 
primary on orbital inclination. 
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readily derivable from spherical trigonometry, if we assume the orbit to be planar. From FIG 
3 we have 

_ J sin 3 (0/- 0 o ) + Vi sin 2 0 Q sin 2 Q { (1 -xos Q { )' , 

tan I = ; \3') 

cos Oq cos 0 { sin 


For 0 C = 0 (equatorial launch) we have 


tan 1 = 


tan Of 
sin (f>f 


(38) 


A more exact equation can be derived for non-planer trajectories by noting that the inclina- 
tion is defined as the angle between the (instantaneous) orbital plane and the equator. A vector 
perpendicular to the orbital plane is 

r x v = r 2 (— 6 -f <j) cos 0 Uq) (39) 

while a unit vector perpendicular to the equatorial plane is 

k = sin 0 u r + cos 6 u e (40) 


Then 


cos / = ■ 


0 cos 0 


r x v 


(£2 + ^2 cos 2 0) Vi 


giving 


tan I 


y/H±iL sin 2 0 cos 2 0 
(f> cos 2 0 


(41) 


Equation (37) will be used in the formulation rather than the more exact form given by equa- 
tion (41). This will be justified as an acceptable simplification later in the report. (see p. 22, 
para. 2.) 


SECTION III. VARIATIONAL FORMULATION 


In this section we shall formulate the necessary conditions to maximize the gross weight 
placed in orbit. The problem may be rephrased to read: determine the optimum steering program 
that will minimize propellant consumption for ascent from a given initial point to a prescribed 
set of end conditions. 

The mass of propellant expended may be written as 

m P = m o ~ m i = m o ~ ( m o - rntf) - mtf (42) 

Equation (42) shows that minimizing the propellant consumption is equivalent to minimizing the 
time required to attain orbit, as is intuitively apparent for the case of a vehicle with constant 
mass flow rate. 


Let us now make the following definitions: define J { (i = 1,...,6) as 

, . T cos 8 cos y MG _ ■> n 

J. = p - <- + = - r {co 2 + a 2 cos 2 0)~0 

1 r m 

2 pco 


. • T sin 8 

J 2 =co — — + 


+ o 2 sin 0 cos 0 = 0 


. T cos S sin y 2po A 

J 3 = a - 2 co a tan 0 + = 0 

mr cos 0 r 


(43) 

(44) 

(45) 
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(46) 


} 4 -r- P =Q 

Js - <o = 0 (47) 

] 6 = <j) - a = 0 (48) 

Let Aj (z = 1, . . . , 6) be time dependent Lagrange multipliers and write 

F = A-i }j (sum on z) (49) 


Also, define 

r = m t t + Vj [MG -r f 3 (d)| + < 7 / cos 2 ft 


+ i'. 


y- V sin 2 ( 0, - ) + l / 2 sin 2 sin 2 (1 - cos ) 

cos $ 0 cos sin <£ f 


(50) 


where the v's are constant Lagrange multipliers, 


Let y g denote any member of the set p, co, o, r, 6, <f) and x s either of the control vari- 
ables y, 8. 

We have the problem of finding the steering program, y ( t ), 5 ( t ), which minimizes the 
propellant expenditure required to ascend from a given point on the lunar surface into a pre speci- 
fied lunar orbit. This is equivalent to the probleid stated in the following theorem which is 
proved in reference 2. 

Given 


« = T + f F (x s , y s , y a , t) dt 


The necessary conditions to minimize k are given by 

d F 

— = 0 (all time values) 
dx a 

^ 0 (all time values) 
dy a dt \d yj 

For any y a *s not fixed at t = t 0 


(51) 


(52) 


(53) 


ar 

dy Sl 


For any y a ’s not fixed at t = t t 


a F 


= 0 


a r /a f 


d y* f \py 


= 0 


(54) 


(55) 
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For the final time point, t = t f (never fixed) 

777 (su "“ s) 

Since the form of equations (43) through (48) gives explicit representation 
terms we find that 

££_- a 

By. ~ * 

and 

dt \dy a J » 

Equation (53) may be rewritten in the form 

<• -(??•» 


SECTION IV. VARIATIONAL EQUATIONS 

The expression for F, defined by equation (49), is given in expanded form by 
. r . T cos 8 cos y MG / 2 , 2 

p - x > L p - K- + “ " rU +a cos fl) J 

+ \ I" <ii - 7 - sin 4 7 + + CT 2 sin 0 cos tf] 

2 [ (m 0 - m t) r r J 

. f . T cos § sin y 0 . ^ % P a l 

+ A, CT - 7 . , r /3 — 2 (i) <7 tan o + — ^ — 1 

3 L (772 n - mt)r cos 0 


r J 


+ A 4 [ r - p ] + A s [ 0 - a> 3 + A 6 

Applying equation (52) to the control variable y 


d F 


dy (ttiq - m t) 

For a constant thrust single stage vehicle 

T 


s [ A > sin * - r^re cos 7J - 0 


( m Q - m t) 


+ 0 


Thus we must have either 


5=0 


(56) 
for the y, 

(57) 

(58) 

(59) 


(60) 

(61) 

(62) 

(63) 
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I I 


or 


Aj sin y 


*3 

r cos d 


cos y 


0 


(64) 


Equation (63) cannot be true at t = t 0 since polar orbits have been ruled out by our 

choice of the coordinate system. At other time points where equation (63) is not fulfilled* we 
have 


tan y 


^3 

X 1 r cos 0 


(65) 


The Euler equation corresponding to the control variable 8 is 

dF T f\ • » *2 * A 3 . o • 

-=r-=r = 7 s— rj A. sin o cos y — — - cos o + ^ sin o sin y 

38 (m n -mt)\_ r r r cos 6 r 


= 0 


A 3 

cos y + ^77- 

' r cos 6 


[A, 

From equation (65) 


sin y J sin 8 = — - cos 8 


sin y a i 


y/ A 2 + (Aj r cos $) 2 ' 


and 


cos y = ± 


Aj r cos 0 


V A 2 + ( X x r cos 0) 2 


Inserting (68) and (69) into (67) and solving for tan 8 gives 

* ± A, cos 6 

tan 8 = ~ 

yj A| + ( A x r cos d) 2 ' 

Since they will be needed later we solve equation (70) for 

. * ± A 2 cos 6 

sin 8 ~ , . . . . =7— 

V (A2 r 2 + A2).COS2 d + AJ 


and 


cos 8 


= ± / A i r 

V "TX2T2 


2 cos 2 0 + A 2 
2+ A2 )cos2 6 + Af 


( 66 ) 


(67) 


( 68 ) 


(69) 


(70) 


(71) 


(72) 


* During the numerical integration of the set of equations that will result, the probability that 
equation (63) will actually be fulfilled is negligible since cos 8 = 0 occurs only for discrete 
points. It is a simple matter to insert a flag into the computer program to indicate that equation 
(63) is satisfied and thus equation (64) is no longer fulfilled. With this understanding (which 
is sloppy mathematics, but standard engineering) we may proceed under the assumption that 
equation (63) is never satisfied. 
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Turning now to the application of equation (59) we consider the variables y m in the order 
p, co, a, r, 0, cf). 


A, — ^ — =* A* — — (A 2 co + A 3 <j) + A 4 — 0 

Op T 


k 

Aa TO 


= A 2 + 2 (A! rco - + A 3 a tan 6) +X S - 


(73) 

(74) 


A, - 


dF 


3- 5^ 


» A 3 + 2 [ ( Aj r cos 0 - A 2 sin 6 ) a cos 6 + A 3 ( <0 tan 6 - + A 6 = 0 


A 4 - * A 4 + A x [— + ( u> 2 + o' 2 c os 2 # )] + ~ 2 - t ^2 40 + ^ 3 ^ 

- 7 ^-T— r g- F A 2 sin 5 + X 3Cosg.in 7 1 0 

U 0 -wi)r 2 L 2 cos <9 J 


(75) 


(76) 


X — = A* — <r 2 [2A, r sin (9 cos # + A, (cos 2 6 - sin 2 0)] 

5 06 si 

+ 2 X 3 a> o sec 2 8 + A 3 tan 9 sec 9 - 0 




(77) 

(78) 


Equation (78) may immediately Be integrated to give 




(79) 


where Cj is some constant of integration. 

SECTION V. SUMMARY OF ANALYTICAL RESULTS 


The equations of motion, (43) through (48) along with the equations of the turn program, 
(65) and (70), and equations (73) through (77) and (79) (which determine the Lagrange multi- 
pliers that are included in the turn program), completely specify the optimal three-dimensional 
steering program of a single-stage constant-thrust vehicle except for initial conditions. While 
the initial conditions specifying position and velocity at t = t 0 as well as desired end conditions 
may easily be stated, the initial values of Ai,...,A 6 are very difficult to determine. Once these 
values are chosen,no degree of freedom is left in the system. The determination of (A^q, ...,(A 6 )^ 
to obtain desired end conditions is discussed below. 

The equation^ developed in the preceding sections may be reduced in number by inserting 
the expressions for p, co, a from equations (16), (17) and (18) into equations (73) through (77) 
and eliminating the control variables ( y , 8) by use of equations (65) and (70).* Equation (79) 

* The sign ambiguity involved in the sine and cosine functions of the control variables will 
be resolved by choosing an upward launch, toward the north, in the direction of the lunar rota- 
tion. With these restrictions the positive sign should be chosen in each case. Further considera- 
tions of sign choice will be given in Section VI. 
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will be used to eliminate X 6 from equation (75). 

Before making the above substitutions, one further step will be taken to eliminate m 0 from 
the equations of motion. 0 

The term for thrust to mass ratio may be written 


T 


rriQ — mt 


(T/m 0 ) 



(80) 


Now the mass flow rate, m , is related to specific impulse, / , by the equation 

. T 

(go)*Is P 


(81) 


where (g 0 ) e is the acceleration of gravity at the Earth’s surface. Substituting equation (81) 
into equation (80) we have 


m 


o 


_T 

— mt 


(T/m 0 ) 

T t_ 
I S p m Q 


(T/m 0 ) 


m o (8 q >9 ! sp 


(82) 


Dividing numerator 
Moon’s surface, (g 0 >) f , 


and denominator of equation (82) by the acceleration of gravity at the 
we have 


T 


m Q — m t 


__ [ T/ttiq ( g 0 ) t ] 

“ J. T t 

( go\ ™ 0 (go)< (go)® 1 „p 


(83) 


The term [ T/m Q (g Q ) t ] can now be recognized as the initial thru st-to -weight of the vehicle 
evaluated at the lunar surface. We shall abbreviate this" term as a. Thus 


m 


o 


r 

- m t 


1 

(go\ 


a 


(gO ^ 3 p 


a (gp)i. (gq)® 1 mp 

(gO )® 1 Bp “ a (gO )< 1 


(84) 


The right-hand member of equation (84) is preferable to the left-hand member since it shows 
direct dependence on the important parameters of initial thrust~to -weight and specific impulse. 
The left-hand member masks this dependence by the inclusion of the lift-off mass. 


Performing the above indicated modifications, we now summarize the results of the pre- 
ceding work with the following set of equations as the final form of equations (13) through (15) 
and (73) through (78). 


r a(go)< (g ojffi Kd 1 
t(g 0 )<B ! b P ~ a (g 0 )< 


+ '■<'<? 2 +0 2 cos 2 0) 
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q[ a(gg>* (gq)* K* ~ [ 
USo^ Z -p ~ a( 8o)< 1 -* 


2 tQ 


— ^> 2 sin 9 cos 9 


( 86 ) 


A 2 cos $ 


r t/(A£ r 2 + A 2 ) cos 2 0 + A^ 1 r 
A; 


a = T a i&k (Moh.hs-3 A i = === =r + 20 0 tan fl - 

^ L (g 0 J® f. P - « <W« *“* rcos ^V^5 r2+X2 2 ;cos20+X 3 

X t - ~ (a 2 0 + a 3 <£) + a 4 =o 

A 2 + 2 (Aj r 9 — - 2 — + A 3 0 tan 0) + A s = 0 

A 3 + 2 [(Aj t cos 0 - A 2 sin 0) <£ cos 9 + X 3 (0 tan 9 - r/r)] + Cj* 0 
X 4 + X 1 [^ + (0 2 + <£ 2 cos* 0)] + -ij- (A, e +A,i) 

■]-° 


2r0 

r 


f (80^9 I aa \ 

(A 2 cos 2 6 + A 3 ) 1 

L(g 0 A& l ap ~ a (8oi 1 J 

„ r 2 cos 6 \/ ( A 2 t 2 + A 2 2 ) cos 2 6 + A \ J 


(87) 

( 88 ) 

(89) 

(90) 


(91) 


A* s — <£ 2 (Ajr sin 2 6+ A 2 cos 2 6) + 2 X^d <f> sec 2 # 


+ 


a (go A (gofeUp . — "I 
- (80)9 Kp ~ a (io>* t -J L 


A 3 tan 6 sec 6 

r yj{ r 2 + X 2 2 )co~s 2 6 + Af 



(92) 


The above equations can b? further reduced in number by^ differentiation of equations (88) 
and (89) followed by insertion of the expressions for A 4 and A s . This alternative form is shown 
in Appendix A. 


SECTION VI. CHOICE OF SIGN CONVENTION 

Equations (68), (69), (71) and (72) indicate that there is a certain amount of freedom in 
choosing the sign convention to be adopted in numerical integration of the equations of motion. 
This choice can be investigated either by consideration of the physical parameters (trigonometric 
functions of the gimbal angles) or the mathematical parameters (Lagrange multipliers) with equi- 
valent results. We choose the former method. 

We may consider all possible cases by giving the functions cos 8, sin 8, sin y, cos y either 
plus or minus signs in all possible combinations, and listing the results in tabular form. For 
reasons that will become apparent below, we also list the product cos 8 cos y. 
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Case Number 

cos 8 

sin S 

sin y 

COS y 

cos 8 cos y 

1 

+ 

+ 

+ 

+ 

+ 

2 

+ 

+ 

+ 

- 

- 

3 

+ 


- 

+ 

+ 

4 

+ 

+ 

- 

- 

- 

5 

+ 

- 

+ 

+ 

4- 

6 

+ 

- 

+ 

- 

- 

7 

+ 

- 

- 

4- 

+ 

8 

+ 

- 

- 

- 

- 

9 

- 

+ 

+ 

+ 

- 

10 

- 

+ 

+ 

- 

+ 

11 

- 

+ 

- 

■h 

- 

12 

- 

+ 

- 

- 

+ 

13 

- 

— 

+ 

+ 

- 

14 

- 

- 

+ 

- 

+ 

15 

- 

- 

- 

+ 

- 

16 

_ 


_ 

— 

+ 


Insertion of the sign conventions 9 through 16 into equations 85 through 92 show that they 
correspond, mathematically, to cases 4, 3, 2, 1, 8, 7, 6, and 5, respectively. Thus we immedi- 
ately reduce the^number of cases to be considered by a factor of two. 

Furthermore, in cases 2, 4, 6, and 8 the sign of the product cos 8 cos y is negative on the 
radial thrust term. If the initial value of A i were chosen negative then the resultant force would 
be in alignment with the positive radius vector for at least a portion of the powered flight. A few 
numerical experiments were conducted to check this possibility, and it was found that the cor- 
responding trajectories did exist. 

An analagous situation is found in cases 5 and 7. For case 5 we can obtain a positive 
inclination by first setting the thrust term on the 6 equation negative (the signs on the thrust 

terms are, of course, arbitrary) and then forcing a negative value of sin 8. Thus, we obtain 

another possible solution. 

Finally, cases 1 and 3 will be considered. Case 1 is the most physically reasonable of 
possible choices. The first of these two cases is equivalent to firing with the planetary rota- 
tion to the north in an outward direction (i.e., in the direction of increasing radius vector). 
Case 3, with tan 8 negative, corresponds exactly, except that we now fire against the plane- 
tary rotation. 

The choice of which of the above cases to be studied will now be discussed. Cases 1 and 
3 are, of course, of the most importance. Case I is studied in detail and case 3 was given less 
consideration. In order to check the possible solutions corresponding to the other remaining 
cases, case 5 was carried out in equivalent detail to case 1 and case 7 in a manner similar to 

case 3. Cases 2, 4, 6, and 8 were given only a cursory check and the results are not pre- 

sented. 

Another point of consideration before closing this section is symmetry under 6 reflection. 
From physical considerations it is fairly obvious that if the initial latitude is Oq and an orbit 
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of inclination 1 achieved from ascent from this point, the equivalent physical results must be 
obtained if the launch latitude is — 6 0 and the resulting inclination is ■*? /. As is well known, 
the fre e-flight (equations possess this symmetry. We may then require that the powered-flight 
equations, including the equations governing the Lagrange multipliers, possess this same sym- 
metry. 

The problem may be approached by replacing 6 by — 6 in equations (85) through ( 92 ) and 
observing the results. This procedure shows that our principle holds only if the multipliers 
Aj, A 3 , and A 4 are symmetric with respect to 6 reflection, while the multipliers A 2 and A s are 
antisymmetric with respect to this reflection. This result might have been predicted since the 
latter two multipliers are those associated with non-planar flight. 

We may now essentially double the amount of data available from our results by the above 
principle. Suppose, for example, that we have numerical data corresponding to the initial condi- 
tions 6 =0, and a final inclination of 5°. We may then obtain identical results by reversing 

the sign on the inclination and the signs on the Lagrange multipliers A 2 and A 5 . 

In the numerical data presented below, the sign convention used in preparation of the data 
will be specified by use of the case numbers given above, i.e., "sign convention for case number 
1 ** "sign convention for case number 3 ," etc. 

One final remark is in order before proceeding to the numerical integration procedure. A 
comment was made above about "forcing” a negative value of sin 8, and the immediate question 
that comes to mind is just how this may be done in practice. 

More generally we may consider the problem of arbitrarily fixing the signs of the trigonometric 
functions of the two thrust orientation angles. We consider the signs of the sines to demon- 
strate the involved principles. Equation ( 68 ) shows that sin y will be positive or negative 
according to the sign of A 3 if the radical is always taken as positive. Likewise, for cos d posi- 
tive and the radical of equation (71) positive, the sign of sin 8 will agree with the sign of A 2 . 

Thus the problem is reduced to fixing the signs of A 2 and A 3 . The procedure for doing this 
is not immediately obvious. It was found, numerically, that A 5 and Cj exerted such strong 
dominance over all other terms in equations (89) and ( 90 ) that A 2 and A 3 always have the opposite 
sign of As and C 1 , respectively. The large values of As that were necessary to give even low 
inclinations, as well as the large value of C 1 that was chosen, ensured that this sign asymmetry 
was maintained throughout the powered flight. 

In summary, then, reversing the sign of the initial value of As reverses the sign on the final 
value of the inclination with no other changes. Reversing the sign of Cj changes the direction 
of launch with respect to the lunar rotation. In the latter case other parameters also change, 
because there does not exist an -east-west symmetry on a rotating sphere. 

VII. NUMERICAL INTEGRATION PROCEDURE 

Since the resulting set of Euler equations, along with the equations of motion, could not be 
solved by analytical methods it was necessary to resort to numerical integration techniques. 

In order to obtain some feel for the results that can be obtained, we note that at the initial 
point the variables r 0 , t Q , 6 0 , d 0 , <£ 0 , (fi 0 , Aj , A 2 , A 3 , A° , A s , Cj , I gp , and a 

(thrust-t o-weight ratio) must be specified. The first six of these ( r , r, 0 , d), <^>,$>) are fixed by 

our choice of launch site. The final two (I and a ) are fixed either by technological feasibility 
\ . . 

or optimization criterion. 
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The X’s are further restricted by two considerations. Once Xj has been chose*}, the lift-off 
angles y Q and S 0 then determine X 2 and X 3 . Furthermore, one of the multipliers acts only as a 
scaling factor since equations (88) through (92) are homogeneous * 

Thus, we are left with only X* , X 4 , and X 5 as arbitrary initial parameters. 

At the end point there are numerous conditions that it would be convenient to specify. Re- 
stricting our discussion to circular orbits,** the first such requirement is that circular energy is 
obtained. Secondly, the requirement that r t = 0 must be invoked to insure circularity (this 
corresponds to a flight path angle of 90 degrees from the vertical in body-fixed coordinates). 
Thirdly, the requirements on the final orbital inclination must be attained. From stability of 
guidance considerations the thrust vector must be aligned with the velocity vector at injection. 
(This puts endpoint requirements on y/ and 8 f .) Finally, the altitude should be a predetermined 
variable. 

We are apparently faced with the problem of mapping three initial conditions into six end 
point conditions. One further degree of freedom is available, however, namely the choice of 
cutoff time. This parameter serves very nicely to determine circular energy. 

The problem may now be stated as: 

Given r Q , r 0 , 6 0 , d Q , cf> Q , cf> Q , y 0 (orX 3 ), 8 0 (or X 2 ), / sp and a, determine 
*1, *5 and h 

such that cutoff energy, zero radial velocity, specified inclination, specified altitude, and final 
alignment of the thrust vector with the velocity vector are obtained.*** By elementary consid- 
erations it is not possible to map four initial conditions into six end point conditions. 

In practice, circular energy, zero radial velocity and desired inclination are always required. 
This leaves one free parameter, and a choice of the most desirable remaining end point to be 
obtained must be made.**** 


It is interesting to note that there is not a one-to-one control correspondence between initial 
and final conditions. That is, one cannot state that Xj controls final radial velocity, X 4 controls 
final gamma, etc. In reality, the situation is much more complex and there is a complicated 
interrelationship between initial conditions and end conditions. The method employed was, first 
of all, to check which initial parameters could iterate which end conditions. It was found that \ 1 
could iterate final altitude, radial velocity, or final gamma. X 4 was good for isolating final radial 


* The multiplier chosen to be arbitrary was C 1 (or X 6 ). As was noted previously, a negative 
C 1 corresponds to firing with the lunar rotation and positive corresponds to firing against the 
lunar rotation. No orbit could be attained with C 1 = 0. The value chosen for C x was -10 5 . 

** The problem of non-circular orbits will be discussed in APPENDIX B. 

*** It is very important to realize, at this point, that the angles <J> and 6 at cutoff ( <J) t , 8 1 ) are 
fixed (though unknown). This fact is very important in connection with the necessary end point 
conditions. 

****From this point on we refer to isolating final gamma (= n / 2 in all cases) as an equivalence 
to isolating the angle between the thrust and velocity vectors. 
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velocity, less adept at final altitude and very poor for final gamma. Only A 5 was used to de- 
termine inclination— and it was used for nothing else. ' Checks were also made to determine the 
applicability of A 3 and Cj as iteration parameters. A 3 was found to isolate final gamma and 
nothing else. Q can be used as an iteration parameter (if another A is frozen) and was found to 
Be acceptable (but not good) for isolating final radial velocity and final gamma, but very poor 
for isolating final altitude. 

The iteration parameters were then chosen and a ’'forced’' one-to-one correspondence was 
used. For example, A x was chosen to isolate a final gamma of 90 degrees, A 4 was chosen to 
force the final radial velocity to be zero, and A 5 was chosen to determine final inclination. Once 
the desired gamma was obtained in a first order isolation, the computer incremented A 4 , recon- 
verged Ai , etc., until both final gamma and radial velocity matched specified end conditions. 
Then the inclination was attached in a third order isolation. Thus, one-to-one correspondence is 
only a surface artiface and the computer, in reality, searches out a single three-dimensional 
initial point which corresponds to the prespecified final point. This method is generally referred 
to as the "cruddy creeper”. Another method of approach is to attempt a complete run with 
guessed initial multipliers and record the end conditions. One of the multipliers is then modified, 
and the end conditions again recorded. The multipliers are then reset to the initial guesses and 
the above procedure is repeated until all multipliers have ■ been changed and the results noted. 
The machine then performs a multi-dimensional interpolation for the initial conditions which 
yield the desired end conditions. This procedure was finally chosen in preference to the “cruddy 
creeper” since its convergence time is between one and two orders of magnitude faster. 

It was found that there also exist initial values of the thrust-to-weight ratio which 
produce maximum values for the mass fraction, if the final altitude is specified, and maxi- 
mum final altitudes, if the final value of gamma is specified to be 90°. Besides the isolation 
of an ''optimum” thrust-to-weight ratio, it was also found that there exist initial values 
of thd 1 lift-off angles, gamma and delta, which maximize mass ratio for the second sort of 
trajectory. We could thus have trajectories which isolate four end conditions and three initial 
conditions. Such trajectories would require a prohibitively long running time, even on the 
fastest available computers. 

The trajectories which maximize mass fraction via choice of initial thrust-to-weight ratio 
are relatively easy to isolate since the maxima are rather flat and the burning times are short. 
The isolation of the thrust-to-weight ratio which maximizes final altitude (for the case of * 
n/2) are extremely difficult to isolate and these trajectories are quite unstable. The isolation 
of initial values of gamma and delta which maximize the mass fraction at cutoff is also difficult 
but is of less practical importance. 

The choice of the time increment used in numerical integration is a very important quantity. 

A previous version of this report (Ref. 3) presented data obtained by using a value of 64 seconds 
for numerical integration. (These data are included here.) This time increment was initially 
chosen by comparison of trajectories run with varying time increments. Unfortunately, the initial 
thrust-to-weight ratio used for the comparative runs was low and false results as to the relative 
accuracies were obtained. The difference in end conditions has since been found to vary strongly 
with high time increments for trajectories which have a high initial thrust-to-weight ratio. This 
is because these trajectories have shorter burning times and thus fewer segments are used during 
the numerical integration than in the low initial thrust-to-weight ratio cases. 

Thus far the problem is simply one of accuracy. The final values of the several variables 
corresponding to the higher thrust-to-weight ratios were found to be enough in error that false 
maxima were introduced into the resultant graphs of mass fraction (ratio of initial mass to final 
mass) as a function of initial thrust-to-weight ratio for those trajectories which align thrust and 
velocity at cutoff. 
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Restricting ourselves to this case, for the moment, trajectories were studied which maxi- 
mized m i / m 0 as a function of T/W 0 for various values of the time increment used in numerical 
integrations. FIG 4 shows the results of this study. It may be seen that, as the time increment 
approaches zero, the optimum initial lunar thrust-to-weight ratio apparently becomes unbounded. 

This result may be predicted in retrospect. As is well known, the optimal burning program 
to maximize mass fraction (exclusive of atmospheric drag) is impulsive. We might thus expect 
that if the vehicle is launched in the direction of the lunar rotation at an arbitrarily high thrust 
level, we would obtain a maximum mass fraction, since gravity losses have vanished. It is not 
immediately apparent that this argument is applicable since the cases treated assume non- 
horizontal lift-off and it is not possible to launch in the direction of lunar rotation. It should be 
remembered, however, that the equations of motion are derived under the assumption of a point 
vehicle, which is equivalent to assuming vanishing moments of inertia. The turning rate of the 
vehicle is, thus, governed by only the momentum of tlje point mass, which is negligible at lift-off. 

The above result is, then, the consequence of the assumption of a point vehicle in any area 
where this approximation is invalid. The question arises as to whether the mass fraction may be 
maximized via a choice of initial thrust-to-weight ratio if the altitude is prespecified rather than 
alignment of thrust and velocity vectors. Such maxima were found to exist. In this case these 
maxima were essentially independent of the time increment used in integration for values of the 
increment of four seconds or less. 

The actual time increment finally chosen was four seconds. This gave accuracy of at least 
five significant figures by comparison with time increments of one second. The numerical data 
given in. a previous version of this report (Ref. 3) were run at a time increment of 64 seconds and 
are also included. The purpose for this inclusion is to provide a large amount of useful data that 
is sufficiently accurate for preliminary design purposes. A duplication of these data at a lower 
time increment was made only for cases which were considered to be of particular importance. 
Furthermore, the entire set of data which deals with ascent to a 15-kilometer orbit was not in- 
cluded in the previous publication. 

The technological range considered was chosen to give a broad brush outline of what values 
may be of interest rather than a detailed study of a given configuration. The specific impulses 
considered were 300, 350, 400, and 450 sec. The thrust-to-weight ratios chosen were 1, 1.1, 
1.3, 1.5, 2, 3, 4, 5, 6, and 7 (lunar reference).* 

The launch site was arbitrarily chosen to be the lunar equator at the lunar prime meridian 
{<f> Q = d o =0) and the final orbital inclination was chosen to be 5 degrees. 

The lift-off angle, y 0 , was varied from 0° to 40° in steps of 10°. The angle 8 was always 
set equal to zero. For the special case of y 0 =0°, the initial thrust-to-weight ratio was con- 
sidered for the range between 1 and 2. 

A few cases are given in the next section to illustrate the procedure for obtaining orbits of 
inclination higher than 5°, as well as the effects of non-equatorial launch. 

Below is a resume of the constants used in the program. These values were the best avail- 
able at inception of this report, but probably can be improved upon due to more recent measure- 
ments of various constants. 

(MG), = 4.899996 X 10 12 m 3 /sec 2 
(r 0 ), = 1.738 x 10 6 m 

(go), = 1-622169 m / sec 2 

(g 0 )© = 9-81 m/ sec 2 

APPENDIX C shows a block diagram of the final form of the computer program. 

* To convert from lunar referenced thrust-to-weight to a corresponding value referenced to the 
Earth’s surface it is necessary only to multiply by 
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g, /g, = -1653587 



Predicted Optimal Initial Thrust-to-Weight Ratio 


I 



FIGURE 4. PREDICTED OPTIMAL INITIAL THRUST-TO-WEIGHT RATIO 
VS TIME INCREMENT USED IN NUMERICAL INTEGRATION; 
THRUST AND VELOCITY VECTORS ALIGNED AT ORBIT. 
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SECTION VIII. NUMERICAL RESULTS 


In this section, we present the results of numerical integrations performed on the set of 
simultaneous equations (85) through (92). These results, contained in Tables 1 through 52 will 
be prefaced by a short introduction on the assumptions made, and an explanation of the purpose 
of each Table. 

A few general remarks can be made about the Tables: 

1. The initial value of the longitude, <^ , is always set equal to zero. Since the assumption 
of rotational symmetry was invoked in the derivation of the equations of motion, this in no way 
further restricts the validity of our solutions. 

2. With the exception of Table 34, the initial value of the latitude, # Q , is always zero. 

3. The initial ■value of the thrust orientation angle § is always zero. 

4. With the exception of Tables 1 through 11, 35, 36, 37, 50, 51, and 52 the final value of the 
thrust orientation angle y is always tt/2. (The final value of 8 was usually within 2°of the 
value of the inclination angle.) 

5. The value of Ciwas always set at — I0 s for inclinations angles in the rqnge 3 tt/2 < / < tt/2 
and + 10 s for inclinations in the range tt/2 < I < 3jt/2 . 

6. The time increment for numerical integration was 4 seconds for Tables 1 through 21 and 64 
seconds for Tables 22 through 52. 

7. The final value of r is always zero and circular energy is always achieved, (r is also zero 
by the conditions of circular energy if y { = tt/2.) 

8. Tables 1 through 37 correspond to the sign convention of case 1 except for the second half 
of Table 33 which corresponds to the sign convention of case 3. Tables 38 through 52 corre- 
spond to the sign convention of case 5 except for the second half of Table 49 which corresponds 
to the sign convention of case 7. It should be noted that analogues of each of Tables 22 through 
37 are given in Tables 38 through 52 except for Table 37. This Table would be identical 
for either sign convention 1 or 5.* 

9. Altitudes are presented in meters, velocities in meter/sec., etc. 

We proceed, now, to a more detailed discussion of the individual tables. Those remarks 
referenced directly to Tables 22 through 37 apply equally well to the analogous Tables 38 
through 52 

Tables 1 through 11 present data for trajectories which ascend to a 15-kilometer orbit under 
the assumption of vertical lift-off. These tables are the result of numerical integration of the 
equations of motion using a time step increment of four seconds. These tables list the final 
value of gamma at orbit and it may be seen that these values are usually near 90°. Those 
trajectories with initial thrust-to-weight ratio of less than four "hump” during the ascent; i.e., 
the final altitude is not the maximum altitude. 

Table 1 contains thrust-to-weight ratios which maximize the percentage of initial mass 
which achieves orbit for each listed specific impulse. 

Tables 12 through 21 always have thrust and velocity alignment at orbit. For this reason 
the altitudes are not constant, but vary over a range of almost 100 kilometers. (It may be noted 
that, for a given specific impulse, there exists an initial thrust-to-weight ratio which produces 
a given altitude. We could thus specify both a final value of gamma and a final altitude, provided 
the initial thrust-to-weight ratio was not fixed. This is reasonable in that another degree of 
freedom has opened up.) 


* See Section VI for a discussion of sign conventions. 
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This set of tables also provides a basis of comparison for later tables which present data 
from a study of the same vehicles with a larger value of the time increment, used in numerical 
integration. These later tables (discussed below) cover a much wider range of initial conditions 
than is presented in the first 21 tables. 

Tables 22 through 25 present data for initial thrust-to-weight ratios of 1, 1.1, 1.3, and 1.5, 
respectively, under the assumption of vertical lift-off. There is a direct correlation between the 
initial thrust-to-weight ratio and the largest lift-off angle that may be used. In the case of verti- 
cal lift-off, the initial thrust-to-weight ratio may, theoretically, be as low as one. On the other 
hand, for an initial lift-off angle of 40° the initial thrust-to-weight ratio was found to have an 
asymptotic limit at about 1.4. For initial thrust-to-weight ratios of less than 2, only vertical 
lift-off was considered because of practical considerations. The reason for inclusion of these 
low values of thrust-torweight ratios for even vertical lift-off is not immediately apparent. More 
will be said about this in Section X, but for the present it is interesting to note that a maximum 
altitude exists in the vicinity of an. initial thrust-to-weight ratio of 1.2. 

Tables 26 through 31 constitute the bulk of the numerical work. Each of these tables has a 
given initial thrust-to-weight ratio, and contains data for varying the specific impulse and initial 
lift-off angle. 

As was noted previously, there is some thrust-to-weight value (for a given specific impulse) 
that produces a maximum altitude. This statement is limited to vertical lift-off, insofar as it 
is investigated here, but is probably more generally applicable. Furthermore, we must bear in 
mind that there is a constraint placed on the final value of y. Under these restrictions, Table 
32 presents initial thrust-to-weight values that produce maximum final altitudes. 

Table 33 is much more specialized than the preceding tables, in that it presents data only 
for the special case of an initial thrust-to-weight ratio of 2, a specific impulse of 300 sec., and 
vertical lift-off. In this table, the vehicle is assumed to lift-off from the equator and ascend to 
circular orbits of various inclinations. It can easily be seen that as we approach the singular 
point of a 90° inclination, the Lagrange multipliers increase at a frightening rate. For this reason 
the inclination values were not studied between 89° and 91° (the time required to isolate even 
these values was quite high). This table is split in half by a dashed line. The data above the 
dotted line (0° < l f < 89°) have a C t value of — = 10 s , and those below ( 91 °< If < 180°) have a 
Ci value of + 10 s (corresponding to sign conventions of cases 1 and 3, respectively). The physi- 
cal parameters corresponding to an inclination of 90° may be found by interpolation. In the event 
that polar orbits are of importance, it would be easier to rewrite the equations of motion with this 
singularity built into an inclination of little importance, or set up the equations in a non-trigono- 
metric form which avoids singularities. 

Table 34 is a short presentation of data obtained for the same vehicle used in Table 33^. In 
this case, the lift-off latitude was varied, and A 5 was always set equal to zero. 

In Table 35, the final value of y is not constrained to be 90°. The vehicle under consider- 
ation is assumed to have an initial thrust-to-weight ratio of 5» and a specific impulse of 450 sec. 
The reference case from which this table was constructed was entry number 16 of Table 29. 
The initial lift-off angle was varied, and A° was used to maintain the same altitude found for the 
case of vertical ascent. 

Table 36 shows the data obtained when the final mass fraction was optimized with respect to 
the initial lift-off angle for the same vehicle used in Table 35. No restriction was placed on 
either the final altitude or final value of gamma. Since two end points were free, one of the 
initial multipliers (in addition to Cj ) was arbitrary. The multiplier chosen was A®.- Several of 
the trajectories in the middle of this table exhibit the phenemonon of ''humping**. Below these 
values, the vehicle "falls** the whole way into orbit, and the "humping** effect vanishes. The 
increase in mass fraction obtained by optimizing the initial value of gamma is about 1.3 per cent. 
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Table 37 optimizes the initial value of gamma to obtain maximum mass fraction in orbit 
under the restriction that the final altitude was constant. The final value of gamma was, of 
course, unconstrainable since A° was used to isolate final altitude. In each of the two cases 
presented, the increase in mass fraction was about 1 per cent. 

As a final point, we comment on the accuracy of using equation (38) (or (37) for non- 
equatorial launch). It was found that for those trajectories achieving an inclination of 5°, that 
equation (38) introduced an error of less than .02°. Even for the case of a final inclination of 
89°, the error introduced by this approximation was only 0.2°. Since the original data was obtain- 
ed via equation (38), it was not felt that this small error was significant enough to justify recom- 
puting all trajectories. 
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TABLE 1. Optimum Initial Lunar Thrust-to-Weight Ratio 
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TABLE 3. Initial Lunar Thrust-to-Weight Ratio 
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TABLE 7. Initial Lunar Thrust-to«Weight Ratio 
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TABLE 9- Initial Lunar Thrust-to-Weight Ratio = 
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TABLE 11. Initial Lunar Thrust-to»Weight Ratio = 
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TABLE 13. Initial Lunar Thrust-to-Weight Ratio = 1.1 



79426 .60199 j 42.16431 1 .0824576 | -8798.1 

















TABLE 15. Initial Lunar Thrust-to-Weight Ratio = 1.5 



39045 .52145 7.99520 .0290282 -8807.6 
42606 .57010 8.61537 .0302629 -8806.4 
45567 .60993 9.13358 .0312711 -8805.6 
48064 .64309 9.57261 .0321097 -8804.9 














TABLE 17. Initial Lunar Thrust-to-Weight Ratio 
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TABLE 19. Initial Lunar Thrust-to-Weight Ratio = 

0 o = 0o = ° ; f,-5° ;y, = */ 2 



6352 .55539 .43458 .0065027 -8804.1 
6843 .60361 .45855 .0067285 -8803.3 
7243 .64283 .47786 .0069076 -8802.7 
7573 .67470 .49373 .0070530 -8802.3 











TABLE 21. Initial Lunar Thrust-to-Weight Ratio = 
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68X94 .54161 148.467 I .2274 -8765.2 

















TABLE 23. Initial Lunar Thrust«to-Weight Ratio = 1.1 
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TABLE 26. Initial Lunar Thrust-to-Weight Ratio = 
Oo = <f> o = 0 ; h = 5° ; y t = n/2 
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TABLE 27. Initial Lunar Thrust-to*Weight Ratio = 
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TABLE 28. Initial Lunar Thrust-to-Weight Ratio 


II 




TABLE 29. Initial Lunar Thrust-to-Weight Ratio = 
^ 0 - 0 0 ~ ^ J 1 i * 5°;y/ = tt/2 
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TABLE 30. Initial Lunar Thrust-to-Weight Ratio = 
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TABLE 31. Initial Lunar Thrust-to-Weight Ratio 
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TABLE 33. Variation of Inclination 
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TABLE 34. Variation of Lift-off Latitude 
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TABLE 38. Initial Lunar Thrust-to-Weight Ratio = 



62700 .45198 62.650 .11215 7761.8 
69591 .49827 71.645 .12393 7632.0 
75518 .53705 80.039 .13488 7516.7 
80671 .57003 87.884 .14509 7413.6 
















TABI^E 40. Initial Lunar Thrust-to-Weight Ratio =1.3 
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400 65636 .58661 23.068 s05450 8298.1 

450 69622 .62015 24.463 .05651 8267.6 
















TABLE 42. Initial Lunar Thrust-to-Weight Ratio = 
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TABLE 43. Initial Lunar Thrust-to-Weight Ratio 
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TABLE 44 . Initial Lunar Thrust>to- Weight Ratio = 
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TABLE 45. Initial Lunar Thrust-to-Weight Ratio = 


i 


i > '<tvO<tnO'cOMr r )( v JNOOoOrof s '0\vOrl<j , vo 

oiHMcno'oncnNNiriNvooiritoO'Oiflin 

ooou»iTicoooOHorsooHnin^ooMinooo 

r^cooo>r-ii^'<ti-«oeMrN.<t-T-iOcnt^<l-t^O-<i- 

OOO'OOMWfl'OHcMCOOlO'-tNOOONOr-lN 


<J'a'Or^NcoNO\vOrltsNh.NinO\inif|iOH 
oowNm<roo<tOtocn<i-ooocno<tHooo 
tMinO'Ovo^a\cooocoHvo(Ncnocoooimr) 
00000>^OC0C0aiOHM0\mOH0MTi^O'-l 
OOOOi-HOOOt— ItHOOOr-lt— IOOO i—It-H 

oooooooooooooooooooo 


i s >oocoLOOffii^nc^fncocno<Nt f iONff)stroco 

<tmin<i-HcoHincn<t-cocoinniooM»N(>N 

ro VD H (Ji (N Nrl^O'j-COHsJO^tOPIM^N f-. 

r^r^ooooor^ooooa^ooooocr.cTvrHcoc»c^Oi-» 


CO 1- INO[»100in<tNHin^ONMO\^HOOiri 

oioo<tO'Ocotnonn<} , <^iNstiNo<j-'Ovon 

NOcomoo'DiTHNvi’^ovoooHcnioo'HrtiflrN 

<j-mLnLnLo<T'0 , iOOOcnco<l'<t<t‘v£)r s ^i^r»r^- 

irilOLnirilOlOiOO^O^O^OvOvOvOvOvOvOvOvOvO 


OOff\OMvOHn'nN^OO'HrNCS [v OQ m VO 
Nr-(OOOvOl s 'HrOHtO^ONO^rOHOO'DOOM 
ooinOioisvj-HvoON^inOst-vDstOsti^^ 
oooooorNvoo^aioooorN(j><^(^oorNOO^ooN 


oooooooooooooooooooo 

oooooLnmmminoooooinLnminm 

nrocnfOcntnnnnM<i , <t<t<t<! , <t4<t<t4 


oooooooooooooooooooo 

H CM m <f H CM CO <t ^-icvlco-ct H CM n 't 





TABLE 46 . Initial Lunar Thrust-to-Weighc Ratio = 
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TABLE 47. Initial Lunar Thrust-to-Weight Ratio = 
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TABLE 49. Variation of Inclination 
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SECTION IX. GRAPHICAL PRESENTATION OF RESULTS 


To present the data given in Section VIII in a more readily useable form, and to aid inter- 
polation, a number of graphs are presented. Section X will deal with specific methods of using 
the graphs (as well as the tables of the preceding section) for sample calculations. 

All of the graphs presented below correspond to the sign convention of case number 1, unless 
otherwise stated (i.e., FIG 60, 61, 62, 63, 64, and 65). 

FIG 5 and 6 ^re a graphical presentation of the most important parameters of Table 1. The 
initial thrust-to-weight ratios shown in FIG 5 are those which maximize mass fraction in orbit. 
The only remaining independent parameter, specific impulse, is shown as the abscissa of these 
graphs. 

FIG 7 through 14 are plotted from Tables 2 through 11. In each case the initial lunar thrust- 
to-weight ratio is chosen as the independent variable since most other quantities depend strongly 
on this quantity. The specific impulse values (which exert less influence) are used as param- 
eters on graphs having multiple curves. 

FIG 15 through 22 are taken from Tables 12 through 21. These graphs are later repeated for 
corresponding data which were calculated using the large numerical integration interval; e.g., 
FIG 15. .“corresponds” to FIG 25 and 26. It may be noted that the mass' fraction at orbit con- 
tinues to increase with an increasing initial thrust-to-weight ratio in FIG 15, but reaches a well- 
defined maximum in FIG 26. FIG 15 is correct; the peak shown in FIG 26 is introduced by the 
large time step used in numerical integration. 

FIG 23 through 33 cover the tabular data for the case of y Q = 0°. FIG 23 is similar to FIG 5, 
but in this case the thrust-to-weight values are chosen to produce maximum altitudes. FIG 24, 
along with FIG 28, 29, 30, 31, 32, and 33 present Lagrange multipliers as functions of specific 
impulse or thrust-to-weight ratios with specific impulse as a parameter. FIG 24 corresponds to 
the same thrust-to-weight assumption as FIG 23. 

FIG 25 and 26 present altitude and mass fraction as functions of initial lunar thrust-to- 
weight ratio with specific impulse as a parameter; they differ only in the range of thrust-to-weight 
ratio that is covered. 

FIG 27 shows a plot of mass fraction vs final altitude with both specific impulse and thrust- 
to-weight ratios as parameters. The data covered are the same as FIG 26. 

FIG 34 through 53 cover the same data as FIG 23 through 33 (but now y Q = 10°, 20°, 30°, 
or 40°) except that in these cases no data are presented for initial lunar thrust-to-weight ratios 
of less than 2. 

FIG 54 to 59 are plots of mass fraction and final altitude as functions of the initial thrust 
orientation angle (y Q ). In these figures, the thrust-to-weight ratio is different for each graph 
(r/W 0 )= 2, 3, 4, 5, 6, or 7). 

FIG 60 through 65 show the’ behavior of r,~ r 0 , m { /m^ A®, A®, and A® (respectively) 1 as 
functions of final inclination for a vehicle with an initial thrust-to-weight ratio of 2, and a spe- 
cific impulse of 300 sec. The sign conventions of cases 1, 3, 5, and 7 (see Section VI) are* 
presented in these graphs. The sign convention influences the final altitudes and mass fractions 
very strongly as can be seen from FIG 50 and 51. Although the magnitudes of A® , A®, and A® 
apparently are unbounded as we approach an inclination of 90®, there appears to be a discon- 
tinuity at infinity. This could be expected from the fact that C ^ has a jump discontinuity across 
polar orbit inclination. 
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FIG 66 through 71 do not correspond to any of the tables of Section VIII, but were prepared 
from a special detailed trajectory printout. These graphs show the time history of various para- 
meters along the trajectory. It is interesting to note that 6 and </> (as well as their first and 
second derivatives) are "parallel** throughout the flight. This ' ‘parallelism** also occurs for the 
gimbal angles y and §. Note that both y and § slightly "overshoot** their final values just prior 
to orbital injection. As can be seen from FIG 71, inclination increases very rapidly during the 
early portion of the flight when the thrust is turning a small velocity vector. 

FIG 72 presents variations in final altitude, mass fraction, and final inclination as a function 
of lift-off latitude. FIG 73 shows A° and A° for the same variation of latitude. The assumptions 
about the vehicle are an initial thrust-to-weight ratio of 2 and a specific impulse of 300 sec. 
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Initial Lunar Thrust-to-Weight Ratio 
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FIGURE 5. FINAL VALUE OF GAMMA, OPTIMAL INITIAL 

LUNAR -TO -WEIGHT RATIO AND MASS FRACTION 
VS SPECIFIC IMPULSE FOR ASCENT TO 15 
KILOMETER ORBIT ( At=4 SEC) 





Initial Lunar Thrust-to-Weight Ratio, T/W 





FIGURE 8. \j VS INITIAL LUNAR THRUST- TO- WEIGHT RATIO FOR ASCENT 
TO 15 KILOMETER ORBIT (At=4 SEC) 
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FIGURE 9. Xj VS INITIAL LUNAR THRUST- TO- WEIGHT 

RATIO FOR ASCENT TO 15 KILOMETER ORBIT 
(At=4 SEC) 
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Initial Lunar Thrust-to- Weight Ratio, T/W q 

FIGURE 10. \ 4 VS INITIAL LUNAR THRUST- TO- WEIGHT RATIO FOR 

ASCENT TO 15 KILOMETER ORBIT (At=4 SEC) 
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Initial Lunar Thrust-to-Weight Ratio, T/W 
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FIGURE 15. MASS FRACTION VS INITIAL LUNAR THRUST- TO- WEIGHT RATIO {T =a, A t=4 SEC) 
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FIG 23. MASS FRACTION, ALTITUDE, AND INITIAL LUNAR THRUST -TO -WEIGHT RATIOS 
VS SPECIFIC IMPULSE FOR MAXIMUM FINAL ALTITUDES {y 0 = 0 ) 
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ALTITUDE, r -r , KILOMETERS 
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FIG 27. MASS FR 
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INITAL LUNAR THRUST-TO WEIGHT RATIO T / Wo 
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FIG 30 . X4 VS INITIAL LUNAR THRUST- TO- WEIGHT RATIO (/o = 0) 
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INITIAL LUNAR THRUST-TO -WEIGHT RATIO, T/w 0 
FIG 32. X c VS INITIAL LUNAR THRUST-TO -WEIGHT RATIO { y = 0° ) 
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FIG 34. FINAL ALTITUDE AND 





MASS FRACTION 
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FIG 35. MASS FRACTION VS FINAL ALTITUDE 0^ = 10°) 
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INITIAL LUNAR THRUST-TO- WEIGHT RATIO, (T/W ) 0 


FIG *2. X 4 VS INITIAL LUNAR TH RUST - TO - WE IGHT RATI 0 < y Q * 20*) 
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INITIAL LUNAR TH RUST- TO -WEIGHT RATIO, IT/W) ( 
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INITIAL LUNAR THRUST- TO - WE I G H T RATIO, (T/W ) 0 
FIG 46. X, VS INITIAL LUNAR TH R UST -TO -WEI GHT RATIO (y Q =ZO°) 


101 


i 





FIG 47. X 4 VS INITIAL LUNAR THRUST-TO- WEIGHT RATIO ( = 30*) 
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INITIAL. LUNAR THRUST-TO - WEIGHT RATIO, 
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FIG 50. MASS FRACTION VS FINAL ALTITUDE (>=40°) 
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INITIAL LUNAR THRUST-TO- WEIGHT RATIO, ( T / W ) Q 

FIG 51. Xj VS INITIAL LUNAR THRUST-TO-WEIGHT RATIO (y 0 = 40°) 
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LIFTOFF ANGLE, Yq . DEGREES 
TION VS LIFTOFF ANGLE FOR (T/W) 
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I NCLIN ATJON, I f , DEGREES 

FIG 63. \4 VS INCLINATION [(VW)o = 2 AND I S p =300 SECONDS] 


118 


I 



119 


to' 


to* 


10 ‘ 



150 180 

INCLINATION , If , DEGREES 

FIG 65 . X 5 VS INCLINATION [(T/W) 0 = 2 AND l #p = 300 SEC.] 
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FIG 67. ANGULAR VELOCITY <£ AND Q VS FU GHT TIME [_{ T/W ) Q = 2 

AND r sp = 300 SECT] 
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FIG 68. ANGULAR ACCELERATIONS <jb AND 9 VS FLIGHT TIME [(T/W) ( 
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FIGURE 77. INITIAL VALUES OF LAGRANGIAN MULTIPLIERS VS 
SPECIFIC IMPULSE FOR INITIAL LUNAR THRUST - 
TO -WEIGHT RATIO OF 4.4796 
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FLIGHT TIME, 

INCLINATION VS FLIGHT TIME [(T/W) 0 = 2 AND I sp = 300 SECONDS 



FIG 72 . FINAL ALTITUDE, MASS FRACTION, AND INCLINATION VS LIFTOFF LATITUDE 

[( T / W )o=2,I S p = 300 SEC.,X 5 =0] 
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SECTION X. EXAMPLES 


In order to illustrate a few of the possible methods of applying the preceding material to 
practical cases, a number of examples will be given. Since this entire report is unclassified, 
the following numbers do not correspond to real vehicles, but the methods of application carry 
over directly. While most cases below are treated graphically, one could use finite difference 
methods to interpolate directly from the tables. 

Example 1. Suppose that we wish to determine the payload that an engine having 15,000 
pounds of thrust and a specific impulse of 440 sec. can place in lunar orbit as a function of lift- 
off weight. Assume the lift-off weights to be studied are specified to be 25,000 pounds and 
30,000 pounds (earth reference). Require that thrust and velocity vectors be aligned at orbit. 

The first step is the construction of Table 53: 


TABLE 53 


'.p 


T 

cr/w Q ) e 


mm 

■ 

.60000 


Wtism 


.54545 


HI 

1 toil 

.50000 


The last column was obtained by dividing the thrust ( T ) by the initial weight (W Q ). The 
next step is to convert the initial thrust-to-weight ratio from earth-reference to lunar-reference. 
This may be accomplished by taking the ratios of the earth value of the acceleration of gravity 
to the lunar value and multiplying the fourth column by this ratio. 

Denoting the Moon by the subscript < and the Earth by the subscript © we have 



9.81 m / sec 2 

1.622169 m/sec 2 


= 6.047459 


Multiplying the earth-referenced thrust-to-weight by 6.0 47459, our table now looks as follows: 


TABLE 54 


{ a P 


T 

(T/W 0 ) @ 

(T/W„ ), 

440 

M 

15,000 

.60000 

3.6285 

440 

ESssfl 

EQjZjl 

.54545 

3.2986 

440 

mm 

m 

.50000 

3. .02 38 


Now we are in a position to use the material presented above. The bulk of the data given 
assumes launch from the lunar equator at the zero meridian into a plane of 5° inclination to the 
equator. A further assumption is called for about the lift-off angle. The data presented covers 
lift-off angles from vertical to 40° from vertical. While 40° lift-off angles give better performance, 
vertical ascent is probably a more reasonable assumption. With this in mind, we choose vertical 
ascent (y 0 *= 0°). 
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Since the thrust-to- weight is specified, we cannot choose an optimum value. FIG 15 shows 
the mass fraction and FIG 17 shows final altitude as functions of the initial lunar thrust-to-weight 
ratios covering the range of interest with parametric values of 300, 350, 400, and 450 under the 
assumption of vertical ascent and alignment of thrust and velocity vectors at orbit. 

Using this figure, we begin by constructing a vertical line from (T/W 0 \ = 3.6285, (T/W 0 ),= 
3.2986, and (T/W 0 ) ( = 3. 0238. The following table is obtained: 


TABLE 55 


(T/W 0 ) t 

hp 

m f/m 0 

n- r 0 

3.6285 

450 

.6648 

18.2 

3.6285 

400 

.6325 

17.4 

3.6285 

350 

.5930 

16.4 

3.6285 

300 

.5449 

15.1 

3.298 6 

450 

.6625 

21.4 

3.2986 

400 

.6320 

20.5 

3.2986 

350 

.5905 

19.3 

3.2986 

300 

.5420 

17.7 

3.0238 

450 

.6600 

25.1 

3.0238 

400 

.6278 

23.9 

3.0238 

350 

.5880 

22.5 

3.0238 

300 

.5390 

20.7 


We may now plot m { / m Q and r { — r 0 vs. / with parametric values of (T/Wq )= 3.6285, 
3.2986, 3.0238. This has been done in FI G 74. 

In order to obtain our final data, we now simply construct a vertical line from the specific 
impulse of interest, namely 440. Doing this we find: 


TABLE 56 


( T/W 0 ) 

I sp 

m { /m 0 

T i ' r 0 

3.6285 

440 

.6590 

18.0 

3.2986 

440 

.6568 

21.2 

3.0238 

440 

.6543 



24.9 


Inserting the data of Table 56 into Table 54, our accumulated data are 


TABLE 57 


Is p 

W 0 

T 

(T/W 0 )® 

(T/w 0 x 

1 

S' 

0 

0 " 

440 

25,000 

15,000 

.60000 

3.6285 

.6590 

18.0 

440 

27,500 

15,000 

.54545 

3.2986 

.6568 

21.2 

440 

30,000 

15,000 

.50000 

3.0238 

.6543 

24.9 
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FIGURE 74. ALTITUDE AND MASS FRACTION VS SPECIFIC IMPULSE 
FOR INITIAL LUNAR THRUST -TO -WEIGHT RATIOS 
OF 3.0238, 3.2986, and 3.6285 
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The payloa4 may now be found (in earth pounds) by multiplying the initial weight, W 0 , by the 
mass fraction placed in orbit m { /m Q . Doing this, we find the following: 


TABLE 58 


'*p 

(sec.) 

^0 

(pounds) 

T 

(pounds) 

(T/W 0 )e 

(Earth) 

(T/W 0 ), 

(Moon) 

inf /mo 

rf~ r 0 
(kilometers) 

Payload 

(pounds) 

440 

25,000 

15,000 

.60000 

3.6286 

.6590 

18.0 

16,480 

440 

27,500. 

15,000 

.54545 

3.298 6 

.6568 

21.2 

18,060 

440 

30,000 

15,000 

.50000 

3.0238 

. 

.6543 

24.9 

19,630 


The last column of Table 58 gives the payload, in pounds, that can be delivered to the 
altitude listed in the next-to-last column. 

Although this completes the solution, there is one other datum that is of interest. The 
characteristic velocity, u, may easily be derived from the mass fraction m f /m Q . This quantity 
is defined as 


“ = So ! ap ln R 

where R is the mass ratio [i.e ., the reciprocal l/(m { /m 0 ) 

a = go l.p In — - (9.81X440) ln 

r 171 f 

Inserting values o( m { /m 0 from Table 58, we may calculate 


1 

( m { /m Q ) 


w ( /jb 0 ] . We may then write 

1 


= 4316.4 ln 

u as follows: 


(m t /m 0 ) 


TABLE 59 


ruf/mo 

1 /(rn t /m Q ) 

In 

\[v=4 

3 m ol 

u 

.6590 

1.517 

.4167 

1799 

.6568 

1.523 

.4207 

1816 

.6543 

1.528 

.4240 

1830 


The final table, including u is 

TABLE 60 


'.p 

(sec) 

Wo 

(pounds) 

T 

(pounds) 

(T/W 0 ) e 

(T/W 0 ) ( 

m i/m 0 

ri - r 0 

(kilometers) 

Payload 

(pounds) 

u 

(meters/sec) 

440 

25,000 

15,000 

.60000 

3.6285 

.6590 

18.0 

16,480 

1799 

440 

27,500 

15,000 

.54545 

3.2986 

.6568 

21.2 

18,060 

1816 

440 

30,000 

15,000 

.50000 

3.0238 

.6543 

24.9 

19,630 

1830 
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Table 60 is reproduced using the approximate data of FIG 26 for purposes of comparison. 
It is shown below as Table 61. 


TABLE 61 


(sec) 

(pounds) 

T 

(pounds) 

( T/W „ ). 

(r/iy, 

m f /m Q 

r / “ r o 
(kilometers) 

Payload 

(pounds) 

u 

(meters/sec) 

440 

wm 

mm 

.60000 

3.6285 

.6605 

17.6 

16,513 

1790 

440 

EE 1 

msm 

.54545 

3.2986 

.6587 

21.1 

18,114 

1802 

440 

EH 

wm 

.50000 

3.0238 

.6554 

24.9 

19,662 

1824 


Example 2. A second case which illustrates another interesting point of the theory is to 
assume a specific impulse of 395 seconds, and a thrust of 20,000 pounds. Let us consider 
lift-off weights of 27,000 pounds and 33 ,000 pounds. Assume that an orbital altitude of 15 
kilometers is specified. 

As before, we construct the following table 


TABLE 62 


I nn 

sp 

^0 

T 

(T/%)fe 

(T/%\ 

395 

27,000 

20,000 

.74074 

4.4796 

395 


20,000 

.66667 

4.0317 

395 

33,000 

20,000 

.60606 

3.6651 


Constructing vertical lines from (T/W \ = 4.4796, 4.0317, and 3.6651 on FIG 7 we read 

the following values: 


TABLE 63 


' W'U 

1 sp 

(seconds) 

m [ /m 0 

4.4796 

450 

.6688 

4.4796 

400 

.6360 

4.4796 

350 

.5958 

4.4796 

300 

.5460 

4.0317 

450 

.6687 

4.0317 

400 

.6360 

4.0317 

350 

.5958 

4.0317 

300 

.5462 

3.6651 

450 

.6684 

3.6651 

400 

.6357 

3.6651 

350 

.5955 

3.6651 

300 

.5460 
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Interestingly enough, all values are so near the optimum thrust-to-weight value~an extremely 
flat region—that it is very difficult to resolve the three mass fraction curves. We thus study the 
thrust-to-weight value of 4.4796. 

FIG 75 shows a plot of m f /m Q vs / for (7*/%),= 4.479 6. From a vertical line at / = 395 
sec, we read 

m f/ m 0 - .6328 

Carrying through this one case (the other two values of thrust-to-weight may be carried 
through by the reader, if desired) we now find 


TABLE 64 


/ 

a p 

(seconds) 

*0 

(pounds) 

T 

(pounds) 

< r/wx 



Payload 

(lbs.) 

u 

(m/sec.) 

395 

27,000 

20,000 

.74074 

4.4796 

.6328 

17,090 

1772 


where u was calculated from 




= (9.81) (395) In 


1 


.6328 

3874.95 In (1.580) = (3874.95) (.4574) = 1772 m/s 


It should be noted that this u t as in example 1, includes all gravity losses, plane change 
maneuvers, etc. 

Example 3. Suppose next that we have the problem of determining the payload capability 
for final orbital altitude of 15 kilometers and characteristic velocity for a moon lift-off using 
/ = 315 sec. Let the engine thrust remain unspecified, for the moment, but consider lift-off 

weights of 16,000 pounds, 18,500 pounds, and 21,000 pounds. 

Since we’ are essentially designing an engine around a specific impulse, we may as well 
choose an optimum thrust-to-weight for each of the above weights. FIG 5 of Section IX gives 
the data needed. This graph shows optimum thrust-to-weight ratios, final value of gamma, and 
mass fractions as a function of specific Impulse for vertical ascent. These values may be read 
directly by constructing a vertical line from a specific impulse value of 315 sec. 

Applying this procedure to FIG 5, we can directly read the following data; 


(T/W Q \ = 4,023 (Optimum) 


(m { /m 0 ) = .5630 

The (optimum) lunar thrust-to-weight ratio may now be converted to an earth-referenced 
value by dividing (rather than multiplying) by the factor (g Q ) e /(<? 0 )« = 6.047459. (The lunar 
thrust-to-weight ratio must always be larger than the earth thrust-to-weight ratio). 

Doing this we find 



4.023 

6.047459 


.6652 
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The engine thrust may now be found by multiplying the initial thrust-to-weight value by the 
lift-off weight. Thus we have 

TABLE 65 


*0 

(pounds) 

(T/W 0 ) e 

T 

(pounds) 

16,000 

.6652 

10,640 

18,500 

.6652 

12,310 

21,000 

.66 52 

13,970 


The payload may now be found by multiplying the mass fraction (the same in each case, 
and equal to .5630) by the lift-off weight. Doing this gives 


TABLE 66 


Wo 

(pounds) 

Payload 

(pounds) 

16,000 

18,500 

21,000 

9,008 

10,416 

11,823 


Since we have one mass fraction for all lift-off weights, we have only one characteristic ve- 
locity which is 


u = g Q l ap In [ l/(m f / m Q )] * (9.8l)(315) In [T/ .5630] = 1775 m/sec. 

Example 4. Consider the problem of delivering a payload to a prespecified altitude with 
thrust and velocity vectors aligned at orbit. Let us assume a specific impulse of 375 sec. and 
an altitude of 30 km. The initial thrust-to-weight ratio is necessarily unspecified since we have 
fixed the final altitude and final value of gamma. 

We begin by construction of a horizontal line at an altitude of 30 km on FIG 16. This line 
then cuts the various altitude curves which correspond to various specific impulses at a definite 
thrust-to-weight ratio. The mass fractions for these thrust-to-weight ratios, and the same speci- 
fic impulse values, may then be read as usual from FIG 15. We obtain the following data: 


TABLE 67 


r f - r 0 (km] 

V sec) 

(T/W 0 1 

m { /m Q 

30.0 

300 

2.37 

.5295 

30.0 

350 

2.52 

.5810 

30.0 

400 

2.63 

.6260 

30.0 

450 

2.71 

.6560 


FIG 76 shows a plot of initial lunar thrust-to-weight ratio and mass fraction vs specific 
impulse for a final altitude of 30 km. 
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The information we desire may now be obtained by constructing a vertical line from a speci- 
fic impulse value of 375 sec., and reading the mass fraction and initial thrust-to-weight ratio 
from the vertical scales. This procedure yields 

(T/W 0 ) = 2.58 
m { /m 0 ss .6042 

This result shows that any vehicle with an initial thrust-to-weight ratio of 2.58, and a 
specific impulse of 375 sec. will place a mass fraction of .6042 into an orbit of 30 km, altitude. 

If the initial engine thrust was specified to be, say, 10,000 pounds, then 

(W 0 ) = = 3875 lunar pounds 

A JO 


The corresponding earth weight is 


(W 0 )© = 23,430 earth pounds 

The payload can now be found, in earth pounds as 

Payload = (W 0 )© (m t / m Q ) = (23,430) (.6042) = 14,160 pounds 

Example 5» Example 3 considered the problem of designing an engine of prespecified 
specific impulse to place a maximum percentage of a given lift-off weight into orbit. Suppose, 
no^, that an engine thrust and specific impulse are given, and we wish to place a maximum 
payload into lunar orbit by variation of the lift-off weight. Suppose, finally, that thrust and 
velocity vectors must be aligned at orbit. 

This problem, at first sight, appears to be equivalent to Example 3. Further consideration 
shows that, a basic distinction exists. The reasoning is as follows: If we fix the specific 

impulse and thrust, and there is a certain lift-off weight that produces the optimum thrust-to- 
weight ratio and maximizes the mass fraction placed in orbit; if we load the vehicle more heavily 

at lift-off, then the final mass fraction decreases, but to find the payload, we multiply the mass 
fraction by the lift-off weight. Although the mass fraction is decreasing, the product of mass 
fraction and lift-off weight (which is increasing) may be either increasing, decreasing, or even 
constant as initial thrust-to-weight ratio varies. 

We propose to investigate this for a specific example. Suppose that the thrust is fixed at 
15,000 pounds, and the, specific impulse is chosen to be 350 sec. (this value can be read di- 
rectly from the tables and will serve for purposes of illustration). 

Let us choose our initial values of the lift-off weight, in lunar pounds, such that the initial 
thrust-to-weight ratios are covered in the tables of Section VIII. The following table is obtained: 


TABLE 68 


/ 

»P 

T 

(pounds) 

ou 

(pounds) 


350 

15,000 

15,000 

1.0 

350 

15,000 

13,636 

1.1 

350 

15,000 

11,538 

1.3 

350 

15,000 

10,000 

1.5 

350 

15,000 

7,500 

2.0 

350 

15,000 

5,000 

3.0 

350 

15,000 

3,750 

4.0 

350 

15,000 

3,000 

5.0 

350 

15,000 

2,500 

6.0 

350 

15,000 

2,143 

7.0 
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The mass fractions may now be read directly from Tables 12 through 21 of Section VIII by 
choosing the entries with / * 350. This yields 


TABLE 69 


(%)< 

(pounds) 

<r/w 0 ) ( 

m // m 0 


15,000 

1.0 

.4729 

7094 

13,636 

1.1 

.4981 

6792 

11,538 

1.3 

.5286 

6099 

10,000 

1.5 

.5465 

5465 

7,500 

2.0 

.5701 

4276 

5,000 

3.0 

.5888 

2944 

3,750 

4.0 

.5967 

2238 

3,000 

5.0 

.6009 

1803 

2,500 

6.0 

.6036 

1509 

2,143 

7.0 

.6054 

1297 


where the last column was found from taking the indicated product. This column shows that the 
payload is continually increasing with decreasing thrust-to-weight ratio, and terminates in an 
end point maximum at (T/W Q ) t «■ 1. 

The implication of the above Table is that a given engine should always be loaded until 
the thrust-to-weight is one to achieve maximum payload. This is a different problem than 
achieving maximum mass fraction.* 

The question of the usefulness of the "optimum’* initial lunar thrust-to-weight ratios now 
arises, and we continue our investigation along slightly different lines. Since there exists no 
optimum for this case, we proceed by choosing an initial thrust-to-weight ratio of 7 for com- 
parative purposes. 

The problem may be approached by considering that the payload was prespecified at 7094 
lunar pounds. Using the above engine, we would then find a necessary initial lunar thrust-to- 
weight ratio of 1. 

But now consider what happens if the above requirement on choice of initial thrust is re- 
laxed, and we are able to choose another engine of arbitrary thrust and the same specific 
impulse. The payload must still be 7094 pounds. 

Table 49 shows that for a specific impulse of 350 sec. the mass fraction corresponding to 
an initial thrust-to-weight ratio of 7 is .6054. The lift-off weight must then be 


ou 



072,/otq) 


7094 

.6054 


11,718 pounds 


T - (T/W Q ) t (W Q ) t = (7) (11,718) - 82,026 pounds 

So far we have merely attained the same payload using two different engines. Now let us 
consider the propellant expenditure used in each case. For the 15,000 pound thrust case, we have 

(% top \ - 0U - <- w t\ - 15 » 000 - 7 » 0 ? 4 “ 7 ’9° 6 P ounds 


* The reader must bear in mind that the earlier remarks about maximization of payload by tra- 
jectory shaping are not modified by this discussion. 
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While for the 82,026 pound thrust case 

(W prop ) 2 = (W 0 \ -:( W t \ = 11,718 - 7,094 - 4,624 pounds 

The difference in propellant expenditures is 

^prop^i - 13 3282 P ounds 

This figure, converted to earth pounds, is 

(W prop)i * 19,248 (earth) pounds 

This difference is the price that must be paid, in this particular example, for a poor choice 
of engine. Since all propellant used for ascent from the lunar surface must first be propelled to 
escape from the earth and braked onto the Moon, the cost is even higher than this example illus- 
trates. 

The above example illustrates quite markedly the difference between various maxima. For 
the first thrust level considered, the payload placed in lunar orbit was a maximum; the second 
thrust level placed the same payload in orbit with a better thrust (referred to mission viewpoint). 
Similar results are obtained for a comparison between the optimal thrust-to-weight ratio in com- 
parison to a thrust-to-weight ratio of unity for the case of ascent to a prespecified altitude. 

It might be objected that the higher thrust engine would weigh a good deal more, and thus 
offset much of our gain. Consider that both engines have a thrust-to-engine-w eight of 25H. 
In practice, we can do much better on the higher thrust engines (aside from clustering), but dis- 
regard this to obtain an upper bound. 

Thus, in the first case we have an (earth referenced) engine' weight of 

W engine = - 600 pounds 


while in the second case 


11/ 82,026 ,- D1 j 

W engine = ^ = 3281 pounds 

Thus, the upper bound of weight gain due to our engine is 2681 pounds while the propellant 
that has been saved is about an order of magnitude greater. 

Example 6. The above cases have been determined from graphical data or read directly from 
the prepared tables. It will often happen that these methods do not yield sufficiently accurate 
results or the required data falls outside the scope of the material presented. 

We shall illustrate the two preceding situations as follows: Suppose that an improved value 
of the final results of Example 2 are required. The initial thrust-to-weight ratio is 4.4796, and 
the specific impulse is 395 sec. The initial problem is to obtain initial values of the various 
Lagrange multipliers which correspond to these initial conditions. 

FIG 9 t 11» a nd 12 show AJ, A£, and A§ (respectively) as functions of initial thrust-to-weight 
ratio. Constructing a vertical line on each of these graphs from a thrust-to-weight ratio of 4.4796, 
the following data are found: 
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TABLE 70 


'.p 

A? 

A°4 


I 

mm 

M 

-8805.7 

-8804.7 

-8803.9 

-8802.3 


FIG 77 shows a cross plot of these data as a function of specific impulse. A vertical line 
from the / value of 395 sec. now yields the values 

A° =1.28 
A° = .01209 
A° =-8803.9 

which can be used along with the values I gp = 395 and T/W 0 = 4.8380 as input data to a com- 
puter (as well as A 2 = A 3 = 0 and C t = — 10 5 ). 

Another problem of interest is to vary the inputs to obtain higher values of the inclination 
than 5°. For instance, we might wish to launch a vehicle into an orbit of 40° inclination. The 
guessing procedure might proceed as follows: Table 33 and the corresponding entry of Table 26 
show that, for the vehicle used. A® increased by 


0 ^ 

<A°,V 


10.44 

7.99 


1.31 


To obtain a 40° inclination orbit we could make a first guess at A° as 


(A° ) 40 ° - (1.31XA? ) 5 ° 


with similar scalings for the other multipliers. 


Similarly, if we wish to vary the lift-off latitude from 0° to 60°, Table 34 gives a scaling 
factor of 


■W >60° _ 15.9?? , n 

(A° ) 0 ° ~ 7.963 " M 

Thus, we could make a first guess at latitude variation of A° as 

) 60 ° =(2)(A?) 0 o 

etc. 


An alternative procedure would be to converge the initial guesses (as was done above) and 
then iterate in small steps for final inclinations; for example, iterate for 10°, 20°, 30°, and 40° 
with the converged initial values for each case fed in an initial data for the following case. 
Similar procedures can be used to extend the ranges of lift-off latitude, thrust-to-weight ratios, 
specific impulses, etc. 
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FLIGHT TIME, t, SECONDS 

FIG 69 . ANGULAR COORDINATES <^AND 9 VS FLIGHT TIME f(T/W) 0 =2 AND- I = 300 SEcTI 


The above guessing game probably appears quite crude to anyone who has not attempted 
an actual isolation of the initial values of the Lagrange multipliers for a case of interest. 
Those who have worked in this area will find it more acceptable. 

SECTION XI. CONCLUSIONS 

The preceding material has covered, of necessity, a rather limited range of problems. This 
restriction is due to the requirement of obtaining all solutions numerically rather than analyti- 
cally. The material presented will, hopefully, aid in initial estimates of the various multipliers 
necessary to obtain other cases of interest. The restriction of constant thrust is contrary to 
most material dealing with variational trajectory shaping, but is more acceptable to present day 
state-of-the-art considerations. 

The necessary end-point conditions were stated symbolically and pursued no further since 
various engineering constraints usually leave no free end point. However, as was demonstrated 
in the numerical results, it is possible to maximize the mass fraction placed into orbit (for the 
case of prespecified orbital altitude) or maximize final altitude (if an angle of attack is specified 
at cutoff) by choice of the initial lunar thrust-to-weight ratio. 

The orbit of maximum mass fraction is of interest from both academic and engineering con- 
siderations. This study was originally undertaken to verify the existence of this orbit. The 
problem is subtler than previously stated (Ref.3) in that a false maximum is predicted if the final 
angle of attack is specified to be zero; the predicted optima, however, depends upon the numerir 
c.al integration step size. On the other hand, prespecification of the altitude that must be 
achieved at orbit predicts a thrust-to-weight ratio that maximizes the mass fraction injected into 
orbit independently of the numerical integration step size; 

The orbit of maximum altitude occurs for low values of the initial thrust-to-weight ratio if 
the angle of attack is specified to be zero at cutoff. These trajectories are quite difficult to 
isolate due to the long burning times and instability with respect to the initial values of the 
multipliers. For thrust-to-weight ratios of less than four (or thereabouts), the altitudes for this 
type of trajectory become prohibitively low. 

The cases which investigate an optimal value of the liftoff angle gamma produce very low 
orbits for a very small increase in mass fraction. Attempts to optimize mass fraction with re- 
spect to the initial value of the angle delta produced extremely unstable trajectories. 

The first two sets of tables are more accurate than the following data since the integration 
stepsize is smaller. It is interesting to note that the later tables are internally consistent. 

A final point with respect to the choice of initial thrust-to-weight ratio is that maximum 
mass fraction does not correspond to maximum payload (see example 5» Section X). The thrust- 
to-weight ratios that are labeled "optimum** refer to the overall mission viewpoint and not to 
the local viewpoint. 

The problem of a choice of sign convention was numerically investigated for four cases; 
two of these four were treated in detail. The most important of these four is the sign convention 
of case number one (see Section VI). This particular sign convention is most important from an 
engineering viewpoint. 
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APPENDIX A 


DERIVATION OF SECOND ORDER EULER EQUATIONS 

The system of five first-order differential equations for Ai > A 2 > A 3 , A 4 , A 5 may be converted 
in a system of two second-order equations (for A j and A 2 ) and one first-order differential equa- 
tion for A 3 by either of two methods. 

The first of these methods involves neglecting the kinematical substitutions for r, 6, <f> and 
writing 


F= Ai [r - 
+ A 2 [§ 

+ A 3 j^ 

The equations 


T cv MG r 2 . Jjl 2 n\ 

— * — ■. cos 5 cos y + — — r \d + d> cos v) 

-mt) ' r 2 ^ J 


(m Q - mt) 


T . * 2 r0 ; 2 . . 

- 7 t ■ ; — sin 8 + + 6 1 sin d cos d\ 

(m n - mt) t r r J 


T • • _ 2 r <h I 

t yr- cos 8 sin y — 2 0 d> tan 6 + 

(m- 7 ht)r cos V r r r J 


for A lt A 2 j A 3 are now given by the second-order Euler equation 


(A-l) 



(A-2) 


(A-3) 


Thus, equation (A 1-2) becomes 





(A-4) 


Note that 


dF 

~cf<jr 


= 0 


Thus, the Euler equation for A 3 becomes 


^3 


d f dF 
dt \ d <f> 


d_ 

dt 


£ - 


= 0 


or 


dF 


(A-5) 


(A -6) 
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A second method of obtaining the second-order equations for A t and A 2 is to differen- 
tiate equations ( 88 ) and ( 89 ) with respect to time, and eliminate A 3 , A 4 , A s , r, 0, <£ from the 
resulting equations by the use of equations (90), (91), (92), (85), ( 86 ) and (87). 


By either of these methods we find the following: 



A 2 6 — A j 


2 MG 
- r 3 


+ 


(0 2 - 3^ 



2A 2 

r 




2 sin 6 cos 6\ 


2 C l( f> 




a (So)< (go)® Kp f A 3 + A 2 cos* 6 

) 9 Kp “ a (s Q h t ]r 2 cos6 LVU 2 ! t 1 + A 2 2 7 cos^ 0 


A* + Ao cos 2 $ 


o'© op 


+ *i 


= 0 


X 2+ 2 (x.rfl- ^-)+X 2 + ^+^(4 sin’ 0-1) -20*] 


(A-7) 


-2A 1 [r(9+2r <£ 2 sin 6 cos 6] + 2 A 3 <£ tan 6 ( 0 tan 6 - r/r) 

A 3 tan 6 sec 6 


-2 C, 4 tan a 4- f -1 - ,, x2 X3 

[(g 0 )© (g 0 ), «J.Tr L VTA^T 


2+ A2 ) COS 2 6> + A 2 


•]-o 


(A-8) 


A 3 + 2 [(A^ cos 0- A 2 sin 0) <£ cos 0 + A 3 (0 tan <9 - r/r)] + C 1= 0 (A-9) 

The set of equations listed above may be used as an alternative set to equations ( 88 ) - (92) 
in determining A lf A 2 , A 3 (the only A*s that appear in the equations of motion). The initial 
value problem is no simpler in this case since we now must guess values of Aj, Ai, A 2 , A 2 ,A 3 ,Ci 
to begin integration. 
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APPENDIX B 


FREE-FUGHT TRANSFER 


Integration of the equations of motion under the assumption of zero thrust will be carried 
out in this Appendix. Although this assumption is contrary to the preceding development, it may 
often be desirable to assume that cutoff occurs at some condition other than a circular orbit. We 
might, for example, bum to parabolic velocity via a circular orbit or burn until an elliptical orbit 
of prespecified parameters is obtained. The following equations will aid in the choice of various 
cutoff conditions for such problems. 

A number of the following equations are, in reality, nothing more than the Kepler equations 
in three dimensions, referenced to the equatorial plane. The standard techniques for integration 
of the two-body equations carry over almost directly. For this reason, no detailed developments 
will be included. 

Under the assumption of zero thrust, the equations of motion may be written: 


f = - — + r {6 2 + d > 2 cos 2 0) 
r 2 r 

(B-l) 

0 — — — — (h 2 sin 6 cos 6 

T 

(B-2) 

0 = 2 0 0 tan Q 2 r< t* 

(B-3) 


Multiplying (B-3) by r 2 cos 2 d an exact differential results and upon integration we find 

r 2 cos 2 6 <f> — (B-4) 

where / 1 (an angular momentum) is a constant of integration. 

Solving equation (B-4) for 0 2 , substituting it into equation (B-2) and multiplying the result 
by 2r 4 d another exact differential results. Upon integration of this expression we find that 


1 6 2 + 


-J! . 

cos 2 6 2 


(B-5) 


where l 2 is the total angular momentum per unit mass. 


Solving equation (B-5) for 6 2 and substituting it in equation (B-l), along with the expression 
for <f> 2 derived above, we find 


M G 

r + 


(B-6) 
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Multiplying this equation by f and integrating yields 


li «c 1 13 P 

rr~ T — — - — — -nr — v + 41 


where E is the total energy per unit mass. 

The final integrations may now be undertaken. From equation (B-4) we can write 


r 2 cos 2 6 d <f> 

T 


= dt 


and from this form the operator — as 


h 


Writing equation (B-2) as 


dt r 2 cos 2 $ d<p 


^ JT t 2 77 (0) ] + sin 6 cos 6 = 0 


and substituting for and ,-^y we 


find 


lx _d_ f h r 2 ddl l\ . . . 

r 4 cos 2 0 d Lr 2 cos 2 0 d <f>\ + r 4 cos 4 0 S1D C ° S 


Noting that 


equation (B-ll) becomes 


which immediately yields 


1 d6 d tan 6 


cos 2 6 d<f> dtf) 


d 2 (tan 0) 
d <fi 2 


+ tan 0=0 


tan 0 = A sin {<f> + <f3 1 ) 


where <£j is an integration constant which determines the position of the vehicle in 
amplitude constant A may be readily evaluated at the cutoff point as 

. tan 6 - 

A = — : — 7 - / c , r\ = tan / 
sm(«56 c + 96 x ) 


The other constant of integration, <j) 1 , may be determined by the use of equation (40). 


(B-7) 

(B-8) 

(B-9) 

(B-10) 

3 

(B-ll) 

(B-12) 

(B-13) 

(B-14) 
orbit. The 

(B-15) 


Solving equhticn^(B-7) for 


dt = 


dr 

YM~G TT 


/ 2E+ ZI 

y r r 


(B-16) 


and equating dt to equation (B-8) we find 


4 


2 E + 


m — To 


r 2 cos 2 6 d </» 


(B-17) 


The 6 term appearing in the last equation may now be eliminated by equation (B-14) yielding 


dr 


__ r*d(f) r 2 drf> 

/op ; 2 m q l\ l ~ l x (1 + tan 2 6) ~ l x [1 + A 2 sin 2 (^6 + x )) 

y r r 2 


Thus 




dr 




(1 i± 

J /j [1 + A 2 sin 2 (<£ + 


&)] 


(B-18) 


Equation (B-18) may now be integrated by setting u = 1/r. Carrying out the integrations yields 


— cos 


-1 


ll/r-MG 


l, 


y/(MG) 2 + 2ElfJ 
Solving for r gives 


2 tan 1 [\/l + A 2 tan'^*^)] 

^ VI + A 2 


1 2 
l J2 _ 


MG + (MG) 2 + 2 E ' cos A 2 tan 1 [\A + ^ tan 




(B -19) 


(B-20) 


Rearranging equation (B-20) into standard form 

H /M G 


r - 


1 + V 1+ Mw cos [ T t \/l + A 2 tan ' 1 [ V 1 + Ai " ““ ^ + & »} ' 


(B-21) 


Two of the orbital parameters, the semi-major axis, a , and eccentricity, e, can be immediately 
recognized from equation (B-21) by requiring the orbit to be a conic section. 




2 Eli 

(B-22) 

(MG) 2 

ll 

£1 

^(B-23) 
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MG 
2 E 


(B-2'4) 


The inclination may still be found by applying equation (40) to the cutoff values of <£ c and 

These parameters will be required to study rendezvous techniques. 

If no orbital inclination change is to occur during free-flight transfer, the above system of 
equations reduce to the more familiar Keplerian description; however, the reference plane must 
be chosen as the same plane in which the vehicle achieves cutoff (u e., circular) conditions. 

Equation (B-21) is far more complicated than necessary. We may begin the simplification 
by noting that, from (B-15) 


L 


l 7 


l x \J l + M 

From equations ( B-4) and (B-5), we have 


l x sec 7 


(B-25) 


4 _ \f d 2 + <f) 2 cos 2 d 


l % sec 1 


(B-26) 


cos 2 6 sec 7 0 

Dif ferientating (B-14) with respect to time and substituting the value of A from (B-15) yields 


x _ tan 7 cos ( <f> + 

~ 1 + tan 2 7 sin2 (cf) + <f> ) 

Eliminating 0 in equation (B-26) via equation (B-27) and the function of 6 by 

1 


(B-27) 


cos 6 ~ 


1 + tan? 7 sin? ( <56 + ^ ) 


(B-28) 


yields 


L 

l x sec 7 


1 


(B-29) 


The resultant trigonometric function appearing in equation (B-21) may now be readily 
simplified. We begin by noting that, for any argument, s. 


cos s = 


1 

V 1 + tan 2 s 


(B-30) 


Thus, replacing s-> tan" 1 s we find 


cos (tan* 1 s ) = 


V 1 + t tan (tan” 1 s) J 2 y/ 1 + s 2 

From this ' 

cos { tan" 1 [sec 7 tan (<f> + <fi x ) ] } = 


sf 1 4 * sec 2 7 tan 2 (<£ + <f>i) 
Replacing sec 2 7 by its expression from equation (B-14) gives 

cos [ tan* 1 [ sec 7 tan ((f) + cf> x ) ] } = cos 6 cos ( (p + <f) x ) 


(B-31) 


(B-32) 


(B-3'3) 
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(B-3-4) 


Assembling the results yields the final form of equation (B-21) as 

a (1 - e 2 ) 

1 + e cos 0 cos (96 +- tf} ) 

Equation (B-7) may be regarded as an expression for r in theory only. Numerical work with 
this equation shows that it is almost useless due to loss of significant figures. For this reason, 
we differentiate equation (B-34) with respect to time to obtain a useable equation for r. Thus 

• a (1 — e 2 ) e [ sin 6 cos ( <f> + 0 ) 6 + cos 6 sin ( <f> + <ft t ) ] 

[ 1 + e cos 0 cos ((£> -f t ^)] 2 


r 2 e [ sin d cos (0 + <£., ) Q + cos 6 sin (0 + <£,) <£] 
= a (1 - e 2 ) 2 


Eliminating 6 and 6 as before we find 


. (r 2 (f> cos 2 0 ) e cos Q sec 2 / sin ( <£ 4 - 

r= g ( 1 — e 2) 


(B-35) 


(B-36) 


Substituting for r 2 cos 2 0 from equation (B-4) and a (1 — e 2 ) from equation (B-23) gives 

. ( l sec I \ f MG\ a . . x 

r = e — J — -j ^ —J cos 0 sec / sj.n ((/> + ^ ) 


= e 



cos 0 sec / sin ( <£ + ) 


(B-37) 


by equation (B-29). A convenient form for this equation may be obtained if we eliminate sin 
(0+ by use of equation (B-14). For / ^ Owe have 


f = e 



sin Q esc / 


(B-38) 


For the case of / = 0 equation (B-37) becomes 


e 



sin + 0 X ) 


(B-39) 


Another datum of importance along the trajectory is time. To obtain an expression for this 
quantity which retains significant figures during numerical manipulations, we write equation 
(B-16) in the form 




rdr 

V 2 E r2+ 2 M G r - l 2 


(B-40) 


Substituting expressions for E and from equations (B-24) and (B-23), respectively, gives 


dt = 



rdr 

yj a 2 e 2 — (a — r) 2 


(B-41) 
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The eccentric anomaly, may be introduced by defining 

a — r =t a e cos & 

From this definition comes 

r = a (1 — e cos & ) 

and 

dr * at e sin (t> d 

Substituting these expressions into equation (B-41) we find 

dt ~ e cos ^ ^ ^ 


(B-42) 


(B-43) 


(B-44) 


(B-45) 




(B-46) 


The eccentric anomaly may now be related to a combination of the angles <^5> and 0 as 
follows: Let us define 


cos 77 = cos 0 cos ( 0 + '<f> ) 
The geometrical significance of 77 is shown in FIG 78. 


(B-47) 


k 



FIG 78. Geometrical Significance of 77 
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Equating the expression for r from equations (B-34) and (B-43) and solving for cos rf results 


From this comes 


and 


cos & — e 

cos 77 = jr- 

' 1 — e cos & 


, (1 - e ) (1 + cos g) 

1 + cos rj = — — — 


1 — e cos g 


. (1 + e) ( 1 — cos g) 

1 — cos 77 - — — — 

1 — e cos g 


Then 


1 — cos 77 ' . 77 / 1 + c ' g 

£_ = • tan -t- — / tan 

1 + cos rj 2 V 1 - e 2 


(B-48) 


(B-49) 


(B-50) 


(B-51) 


The preceding equations of this Appendix will be used to illustrate the method of attack that 
may be used to transfer from a circular orbit to another conic section. 

At any point along the powered trajectory that occurs after passing through circular orbit, 
the values of r, f, 0, 0, <£, <f> are known by numerical integration of equations (85) through (92). 
Thus, the function represented by equation (B-7). 


E = ~ [ r2 + r 2(0 2 + <£ 2 cos 2 0)] - 


MG 


(B-52) 


serves to determine, E. The eccentricity can now be specified from equations (B-4), (B-5), and 
(B-22) as 


-J 1 


(MG) 2 

and the semi-major axis as in equation (B-24) 


1 + ( 6 2 4- j> 2 cos 2 0) 


(B-53) 


MG 

2E 


At each point the apogee altitude may be found from 

r a ~ a (1 + e) (B-54) 

and the perigee altitude from 

r p = a (1 -e) (B-55) 

If a “pseudo-Hohmann” transfer from circular orbit to a higher circular orbit is desired, the 
value of r a will probably be prespecified (as well as the requirement that r p > r). Several 
options are possible, such as requiring that the final altitude of the final circular orbit be equal 
to r a . Or, one might require an elliptical orbit of period equal to that of the circular orbit. 

In any case, the powered flight continues until such time as the specified orbital parameters 
are attained. For a direct parabolic escape one would merely require that cutoff occurs at 
E = 0 (r a * 00 ). 
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APPENDIX C 


COMPUTER FLOW DIAGRAMS 
GENERAL DATA FLOW 
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INTEGRATION BLOCK 

n = Typical Time Step and m = 1, 2, 3, 4, 5 
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Compute 


r 3 = r o + K A t 

r 3 = r 0 + r 2 A t 

Q 3 ~ ^0 + ^2 ^ ^ 

<9 3 = ^0 + ^2 A / 
= ^0 + ^2 A 1 
<k 3 = 00 + ^2 A/ 



To Control Block 
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CONTROL BLOCK 



To Integration 
Block 



To Start Search 
Block 


To Input Block 


To Input Block 


To Start Search Block 



MR BLOCK 



George C. Marshall Space Flight Center 
National Aeronautics and Space Administration 
Huntsville, Alabama, 1963 
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